BROKEN RAY TOMOGRAPHY IN THE DISK 
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Abstract. Given a bounded C'^ domain Q C IR" and a nonempty subset E of 
its boundary (set of tomography), we consider broken rays which start and end 
at points of E. We ask; If the integrals of a function over all such broken rays 
are known, can the function be reconstructed? We give positive answers when 
n is a ball and the unknown function is required to be uniformly quasianalytic 
in the angular variable and the set of tomography is open. We also analyze 
the situation when the set of tomography is a singleton. 



1. Introduction 

Consider a ray of light travelling to the right on the real axis. If the intensity of 
the light at a point x is I{x) and the material on the real axis has a non-constant 
attenuation coefficient / > 0, the differential equation 

(1) I'ix) = -f{x)Iix) 

is satisfied. Given the initial intensity /(O) and the intensity at a distance x, we 
may easily find (under suitable regularity assumptions) that 

(2) ,o, (11)=^' 

That is, by measuring the ratio of initial and final intensity of a beam of light we in 
effect measure the integral of the attenuation coefficient over its trajectory. This is 
true even when the trajectory is more complicated than a segment of the real line. 

The fundamental problem of (scalar) ray tomography is to deduce knowledge of 
the attenuation coefficient / from its integrals over a suitable set of trajectories. Dif- 
ferent materials (and different wavelengths) have different attenuation coefficients, 
so information about / can be translated into understanding of the structure of the 
object being imaged. 

Consider a bounded domain 17 in a Euclidean space and a continuous atten- 
uation coefficient / : — M. A piecewise linear path in ^1 with reflections at the 
boundary as dictated by geometric optics (the angle of incidence equals the angle of 
reflection) is called a broken ray. We have a device E, which is a nonempty subset 
of the boundary dVl, from which we may emit light and measure the intensity of 
incoming light. The rest of the boundary acts as a reflector and we consider those 
broken rays which start and end at E, which we call the set of tomography. The 
corresponding set of integrals of / over such broken rays is called the broken ray 
transform of f with the set of tomography E. We will implicitly assume all broken 
rays to start and end at the set of tomography. 
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We ask the following question: If the integral of / is known over all broken 
rays with both endpoints in E, can the function / be fully reconstructed? How 
does the result depend on the domain Q, the set of tomography E, and regularity 
assumptions on the unknown function /? 

In the case E — dQ broken rays are simply intersections of lines with the domain 
fl. By extending the unknown function / by zero to M" \ fi. the problem reduces to 
the classical X-ray transform. This transform is well understood: The answer to the 
uniqueness question is affirmative, and also regularity and feasible reconstruction 
algorithms are known. For these results we refer to the textbooks by Helgason [S] 
and Natterer [9J and the original works of Radon [10 and Cormack [2j. 

Reflections from a boundary of a domain naturally arise in connection with 
billiards. For details we refer to the textbooks by Tabachnikov JTD and Chernov 
and Markarian [IJ. Reflections in connection to the X-ray transform as considered 
here, however, appear to have been studied less. 

Florescu, Markel, and Schotland [1^ have studied the broken ray Radon transform 
in connection with the problem of recovering absorption and scattering coefficients 
in single scattering tomography. Eskin [3] has studied the reconstruction of elec- 
tromagnetic potential in the presence of convex reflecting obstacles. Lozev [8 has 
also studied the numerical aspect of scalar tomography with similar obstacles. The 
question asked above has also turned out to be related to the Calderon problem 
with partial data. |7] 

In this article we consider the special case where is the unit disk £) C M^, and 
give the following answers: 

Theorem 1. /// : Z) — > M is continuous and the set of tomography is a singleton, 
then integrals over broken rays uniquely determine the integral of f over any circle 
centered at the origin. 

Corollary 2. Let B" denote the unit ball m M". If n > 2, f : i?" M. is 

continuous, and the set of tomography is a singleton, then integrals over broken 
rays uniquely determine the integral of f over any circle centered at the origin with 
the singleton in the circle 's plane. 

Full reconstruction is guaranteed for functions that are uniformly quasianalytic 
in the angular variable in the following theorem and its corollary. Uniformly quasi- 
analytic functions are continuous in the radial variable and quasianalytic]^ in the 
angular variable with a uniformity condition given in Deflnition |14[ There is, in par- 
ticular, great freedom in the radial dependence of the unknown function. Assuming 
quasianalyticity in both variables makes the problem much easier, as demonstrated 
in Corollary |6] 

Theorem 3. If f : D ^ R is uniformly quasianalytic in the angular variable in 
the sense of Definitional^ and its angular derivatives of all orders satisfy the Dini- 
Lipschitz condition and the set of tomography is open, then integrals over broken 
rays uniquely determine f everywhere in the disk. 

Corollary 4. Theorem [5| also holds when D is replaced by the unit ball in any 
Euclidean space of dimension two or higher. 



For example, real-analytic functions are quasianalytic. 
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We need not assume any regularity of / at the origin to guarantee reconstruction 
in the punctured disk D \ {0} in Theorems [l] and [sj If we assume continuity at the 
origin, Theorem [l] imphes that /(O) can be reconstructed. 

One class of functions satisfying the regularity assumptions of Theorem |3] are 
functions of the form (in polar coordinates) 

K 

(3) f{r,d)=ao{r)+Y,{Mr)cos{M)+bk{r)sm{M)), 

fc=i 

where K (iN and the functions ak,bk satisfy the Dini-Lipschitz condition. 

We also show in Lemma [TtI that if Theorem |3] would hold for all functions which 
are C°° in the angular variable, then it would also hold for all continuous functions. 
The question whether Theorem |3] indeed holds for all continuous function remains 
an interesting open question. 

We begin by making some general remarks about the nature of the problem and 
presenting the formalism for the broken ray transform in the disk (Section [2| . In 
Section [23] we prove Theorem [T] and Corollary [2] singleton tomography is discussed 
in more detail in Section [sj After developing the necessary tools in Sections [3] and [4j 
we present the proof of Theorem [3] and Corollary |4] in Section [6] 

1.1. Acknowledgements. The author is partly supported by the Academy of Fin- 
land. The author wishes to thank Mikko Salo for many discussions and insightful 
comments regarding this work. 

2. Preliminaries 

2.1. Some initial remarks. Let us first consider the general problem with a 
bounded Euclidean domain fl. The problem is linear in the function /, so 
we may reformulate the uniqueness question of the introduction as follows: If the 
integral of / over all broken rays vanishes, is / identically zero? 

Linearity also causes the problem for a complex-valued / to reduce to the real- 
valued case: the problems of finding the real and imaginary parts of / decouple. 

Decreasing the size of E makes the problem more difficult by decreasing the 
number of broken rays available for reconstruction. If it is convenient to consider 
some trajectories that have reflections also on E, we may do so. Summing the 
integrals over each part of the trajectory starting and ending at E, we may construct 
the integral over a trajectory with reflections on E. Thus the only requirement for 
the broken rays of interest is that their endpoints lie in the set of tomography E. 

The assumption of regularity of dfl is necessary to make sense of the reflec- 
tions everywhere. A natural assumption on the unknown function / : 51 — !■ M is 
continuity. Corollary [6] demonstrates that too strong regularity assumptions make 
the problem very easy. 

Convexity of Q is not required in the following proposition, but we keep the 
assumption to make the proof simpler. 

Proposition 5. Let ft C M", n>2, be a bounded convex domain. // / : — >■ E is 
continuous and the integral of f over all broken rays vanishes, then f vanishes in the 
interior of the convex hull of each connected component of the set of tomography E. 

Proof. We may assume that E is connected by replacing it by any of its connected 
components and prove that / vanishes in the interior of the convex hull co E of E. 
We extend / by zero to M". 
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Let £ > and denote by the function / convolved with the standard moUifier 
with support in B(0,e). The integral of / vanishes over any line that does not meet 
U = (i\coE, whence the integral of vanishes over any line that does not meet 
the convex set B{U,e) = {x e M" : d{x, U < e)}. Since e C^(M"), the function 
/e vanishes outside B{U,e) by Helgason's support theorem Theorem 2.6]. 

Take any point y G intcoi?. For all e < d{y, U) we know that /e(y) = 0. Since 
/ is continuous at y, f^iy) converges to /(y) as e — >■ 0. □ 

Corollary 6. Suppose C M", n > 2, is a hounded strictly convex domain and 
the set of tomography E C d^l is open. If f is quasianalytic in VL and continuous 
in Cl and the integral of f over all broken rays vanishes, then f is identically zero. 

Proof. We may again assume that E is connected. Since E is open and fl is strictly 
convex, intcoi? ^ 0. By Proposition [s] / vanishes in this nonempty open subset of 
il, from which the claim follows by quasianalyticity of /. □ 

Note that in the preceeding proposition and corollary only the non-reflecting 
trajectories were used, although their union is not even dense in D,. 

2.2. Formalism for the broken ray transform in the disk. We use polar 
coordinates (r, -d) in the unit disk Z) C and identify points on the boundary 
Si = dD with the corresponding angle. The notation E = {0} means that E is one 
point at angle zero. 

A broken ray 7 is characterized by the number of line segments n-y (there are 
Ury — 1 reflections), the initial point L-y e E, the final point G E, the angle a-y 
at which a line segment appears when seen from the origin, and an integer winding 
number m-y such that the trajectory condition 

(4) n^ay = — + 2'Kmy 

is satisfied. These five parameters (n^ , m-y , a-y , t-,, , K-y ) uniquely (and redundantly) 
define the broken ray 7. We require that riy > Q but the other parameters may 
be negative. We denote by Te the set of all broken rays from E to E, allowing 
intermediate refiections on E. We denote by = cos(a^/2) and dj = 2 |sin(Q!-y/2)| 
the distance from the origin to the trajectory and the length of each individual line 
segment. We also write 

(5) = {7 e : > N} 

to denote the set of rays which have sufficiently many reflections. 

A trajectory 7 is a (piecewise linear) curve, and we denote its trace by tr(7). 

Let C = C{D; M.) be the space of continous real- valued mappings from the closed 
unit disk, C^^ C C the subspace of functions which satisfy the Dini-Lipschitz 
condition, and _B(r£;,IR) the space of bounded real-valued mappings from F^;. The 
mappings in i3(F£;,IR) need not be linear, continuous, or measurable (the set F^; is 
not equipped with any such structure), but merely bounded. We equip all of these 
spaces with the supremum norm. The Dini-Lipschitz condition states that the 
modulus of continuity of a function / satisfy uif{S) — o(log~^ ^~^) for small S. 
Holder continuous functions, for example, satisfy the Dini-Lipschitz condition. 
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Definition 7. For / G C, we define the broken ray transform of f as the function 
gf -.Te^M. such that 

(6) Gfii) = ^ fdn' = jI^ fh{t))dt, 

when 7 : [0, — ^ D is a broken ray with unit speed. The map Q : C ^ i3(r£;,R) is 
the broken ray transform. We define the unnormalized broken ray transform Q by 

(7) Qfil) = I fdn\ 



Elementary calculations show that Q is linear and continuous. Using a constant 
function as an example, one may show that moreover \\Q\\ = 1. The unnormalized 
version, however, is not continuous for generic E: the unnormalized broken ray 
transform of a constant function is not bounded, as can be seen by increasing the 
length of the trajectory indefinitely. 

2.3. Proof of Theorem 1. We are now ready to prove Theorem [Tl A Fourier 
analytic proof for the special case when / g C^^ is given in Section^ where the 
result is also discussed in more detail. 

Proof of Theorem^ We denote circles centered at the origin by Sr = S{0,r). We 
denote the average of / over such circles by aQ{r) — fdH^. It suffices to consider 
those 7 that have gcd(n^,m^) = 1. 

Denoting f^{r) = {I — ^)a^ + s arccos(z..y/r), we have 

^/(7) = E / U{r,vtir)) + fir,V7{rm^^^dr 

(8) ' - 

yeS,.ntr(7) 

where we have defined the function pz : (2, 1] — > K for any z €E (0, 1) such that 

(9) p.{r) = {{l-z^){l-z^/r^)y''\ 
Thus 

(10) Gf{i)^^Gfh)^ f'pz^r)-^ /(y)d^- 

The function has a clear geometric interpretation: it is the probability distri- 
bution for the distance from the origin on the trajectory 7. Therefore we expect 
that 5/(7) approaches g{z^), 

(11) g{z) = / p^(r)ao(r)dr, 

J z 

when n.y — 7> 00 in a suitable sense. 

Let us fix some 7 G r|,. The set Sr H tr(7) has 2n^ elements for almost every 
r £ {z^,l). For any such r we may write 

(12) Srntr{^)^{xr,...,x^,n 
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and assume that the points xj'^ are numbered clockwise on Sr- We spHt Sr to 
2n-y pairwise disjoint arcs such that each arc corresponds to the angle 7r/n^ 
and xj'^ £ li'^- This can be done since the points xj'^ are pairwise uniformly 
distributed on Sr in the sense that a relabeling x]'^ i— >■ x]^2 corresponds to a 
rotation by an angle 2n/n^. The length of each arc I^'^ is H^{I^''^) = nr/n^, so 
the diameter is diam(/7''^) < nr/n-y. 

Let e > 0. Since / is continuous in the compact set 13, it is uniformly continuous 
and there is (5 > such that — f{y)\ < e whenever ja: — 2/| < S. If we choose 

7 so that > tt/S, then nr/n^ < S for all r e (z^, 1), whence for almost every 
r G {zj, 1) we find 



(13) 

which yields 

(14) 

and 
(15) 

since pz{r)dr = 1. 



i 



fdu' ~ fixD 



2n^ 



E 

yeSr-ntiif) 



fiy) 



< £ 



1^/(7) - 5(^7)1 <£> 



For any r G [0, 1] we can find a sequence (7^)^^ such that 7; € T^^ and 



as i — >■ cx). By the previous estimate 

(16) lim g/(7, 



9{r) 



independently of the choice of the sequence. We may thus recover g from Q f . Since 



(17) 



2^/\ 



"7' 



we may recover the (zeroth) Abel transform of ag from g. The Abel transform ^0 
is injective and and has an explicit inversion formula (Lemma 10 below), allowing 
us to finally reconstruct the function oq. □ 

Proof of Corollary^ The geometry of the ball _B" is such that any broken ray 7 
lays in the intersection of a two dimensional subspace of M" and the ball. One may 
pick any plane containing the singleton and the origin and reduce the problem to 
singleton tomography in . The result then follows trivially from Theorem [l] □ 

3. Integral transforms 

We write our function / e as a Fourier series in the angular variable. The 
following lemma describes the sense in which the series represents the function. The 
lemma is true also for the complex Fourier series. 

Lemma 8. The Fourier series of any function f £ converges to f uniformly in 
the disk: There are continuous functions '■ [0, 1] — ^ M, k >0, and bk ■ [0, 1] — >■ M, 
fc > 1, such that the sequence of functions {fK)K, 



(18) 



/K(r,i9) =ao(0 



K 
k=l 



afe(r) cos(fct?) + bk{r) sin(fcz9)), 
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converges uniformly to f in the disk as K ^ oo. These functions ak and bk admit 
a common (Dini-Lipschitz) modulus of continuity and satisfy aj.(0) — fefe(O) = 
when k > I. 



Proof. We set for every r G [0,1] 



2 r2TT 



(19) ao(r) = — / f{r,{))M 







and for /c > 1 



2tt 



(20) afc(r) = - / f{r,d)cos{kd)di} 

vr 



and 



(21) bkir) = - f{r,^)sm{k^)M 



2Tr 



These functions manifestly have the property that 0^(0) — 6fc(0) — when fc > 1. 

The functions /(r, •), r g (0, 1], admit a common modulus of continuity which 
satisfies the Dini-Lipschitz condition. (Clearly the functions Uk and bk admit essen- 
tially the same modulus of continuity.) The Fourier series of a uniformly continuous 
function satisfying the Dini-Lipschitz condition converges uniformly and the rate 
of convergence depends only on the modulus of continuity |12[ Theorem 10.3 of 
Chapter 2], whence each fK^f uniformly. □ 

Let / e . Using the notations of the above proof, fx f uniformly and 
Q is continuous in the corresponding topology, whence QfK Qf uniformly. This 
allows us to study the Fourier components and bk of / one by one; this will be 
done in Lemma [T2| and its corollary. 

What we will be able to construct from the broken ray transform Q f is some inte- 
gral transforms of its Fourier components Ok and bk as defined in Lemma [8] Since 
the case k = Q reduces to the classical Abel transform, we call these transforms 
generalized Abel transforms. These transforms appear implicitly in Cormack's ap- 
proach [2J to the Radon transform in the plane. 

Definition 9. For fc = 1, 2, . . . we define the kth generalized Abel transform as a 
map Ak ■■ Co([0, 1]) ^ C([0, 1]) such that 

(22) Akf{z) - 2 f Tk{zly)—M^^==dy, 

J z V 1 " y^ly) 

where Tk is the fcth Chebyshev polynomial and Co([0, 1]) = {/ G C([0, 1]) : /(O) = 
0}. 

For fc = we define the zeroth generalized Abel transform (the Abel transform 
up to a change of variable) as a map '■ C{[Q^ 1]) ^([0, 1]) such that 

(23) Aof{z) - 2 f dy. 

Since To = 1, Eq. (22 1 is correct also in the case fc = 0. 

We summarise the main properties of these transforms. Some of the properties 
are based on those of the Radon transform in the plane, which are covered e.g. 
in [5]. 
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Lemma 10. The generalized Abel transforms defined in Definition^have the fol- 
lowing properties with any fc > 0." 

(1) Ak indeed maps Co([0, 1]) fC([0, 1]) for k ^ 0) into C([0, 1]) as defined. 

(2) Ak is linear. 

(3) The transform is continuous and \\Ak\\ < 2 when C([0, 1]) is equipped with 
the supremum norm^ 

(4) Ak is injective. 

(5) We have the explicit inversion formula 

Proof. (1) For / G C([0, 1]) let g{r,i^) = /(r) cos(fcz9). Also, let 7^ be the Radon 
transform parametrized so that 

(25) ng{r,^)^ f gdH\ 



where X^.i? = {x G M.'^ : xi cos-d + X2 s'md — r}. 
A simple calculation yields 

(26) 7^5(r,^9) = A/(Ocos(M). 

Since /(O) = if fc > 0, .g is continuous, and so is TZg. Therefore the above equation 
shows that Akf is continuous. 

(2) Trivial. 

(3) The result follows from the estimate |Tfe(a;)| < 1 and the observation 

(4) The Radon transform is injective on continuous functions, so the observation 
made in part (1) of the proof shows that Ak is injective. We also have the inversion 
formula 

(28) f{r) = n-^[{s, d) ^ cos{k^)Akf{s)]{r, 0) 

using the inverse Radon transform. 

(5) See Eqs. (10) and (18) in [2 . The treatment therein only requires that / is 
continuous. □ 

4. Tools for open set tomography 

4.1. Explicit form for the broken ray transform. Using the Fourier series rep- 
resentation given in Lemma[8]and the generalized Abel transforms we can explicitly 
calculate the broken ray transform in terms of the Fourier components. Lemma [T2| 
gives the basic result for each Fourier component, and its corollary lists the relevant 
special cases we will need. 

To simplify further notation, we define the coefficient 



sin(fc(K-y — L^)/2) 



when ka-y ^ 27rZ 



(29) Sk{l) = <! sin(fcaT,/2) 

.^(-l)("^+i)^+'''"^ when ka^ e 2tiZ 



^Using Holder's inequality it is not difficult to show that A^k '■ ^''([0, 1]) — > ^'([0. 1]) is contin- 
uous for any p £ (2, oo] and q £ [1, oo]. 
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for all fc = 0, 1, . . . and 7 G F^;. This coefficient has the following properties. 
Lemma 11. For any k and 7 we have: 

(1) \Skil)\<n^. 

(2) The sum identities 



cos(/c(Z — 1/2)0;^) — Skij) cos{k{Kj — i^)/2) and 

(30) 

^sin(A;(Z - l/2)a^) = £'^(7) sin(fc(K^ - t^)/2) 

(3) If {k^ — i-f)/'^ or a^/n is irrational, Skij) ^ 0. 
Proof. (1) By Eq. Q 

(31) sin(fc(K^ — L-y)/2) = sm{knyaj/2 — km-yir) ~ (— 1)'^'™^ sin(fcn-ya-y/2), 
so it suffices to show that the function 

(32) /(X) = ^ 

satisfies |/(a;)| < n for all x e M \ ttZ when n e N and n > 2. 

If we define firmi) = n(— 1)'^"^^^™ for all m S Z, the function / becomes real 
analytic on the real line. In particular, we have /(O) — n. 

By symmetry, it is enough to show that |/(x)| < n for all x G [0,7r/2]. To this 
end we wish to show that /(O) = n, /' ^ on (0, 7r/2n), and |/(a;)| < n when 
X e [7r/2n,7r/2]. 

Indeed, if f'{x) = 0, then ntan(a;) = tan(nx), but this equation has no solutions 
on (0, 7r/2n). And since the sine function is concave on [0, 7r/2], we have sin(7r/2n) > 
1/n, whence l/sin(a;) < n for all x € [7r/2n, 7r/2]. This yields the desired estimate. 

(2) Taking real and imaginary parts of the geometric sum 

{in _ pikP] 
^ whene''^"^! 
ne*'^"/^ when e''=" = 1, 

when a and /3 satisfy na = /3 + 27rm for an integer m, immediately gives the result. 

(3) By Eq. Q (k-^ — i'y)/7r is irrational if and only if a-^jix is. The conclusion 
follows thus immediately from the definition of Sk{'j)- □ 

Lemma 12. Suppose 7 e is such that — 0. If f{r,d) ~ ak{r) cos{ki3) + 
bk{r) sin(A;'i9), where fc = 0, 1, . . . and the functions at and bk are continuous, then 

(34) g/(7) = Sk{j){cos{kK^/2)Akak{z^) + sm(kK-,/2)Akbk{zy)). 

Proof. Let us only consider the case bk = 0; the calculations for the bk term are 
essentially the same. 

Fix some 7 G F^ with l,^ = 0. For any radius r e [0, 1) and direction ip e dD 
we define the corresponding secant Hr^ip = {y £ D : ip ■ y =^ r}. The trajectory 7 
is composed of n-y such secants. Each of these secants corresponds to the radius 
r = z^. We identify the directions tp with angles and label them by ipi = (/ — l/2)ck.y, 
where I — 1,2, ... , n^. 



10 



JOONAS ILMAVIRTA 



After a change of the variable, we find 

f r.nn r^^^ /cos(a../2)\ , ,,cos(a^/2), 

Jh,^,^^ \ COsif) J C0S2((^) 

Since 

(36) cos(fc((^; + (fi)) = cos{kLpi) cos(fc(p) — sin(fc(^/) sin(fc(^) 

and the second term vanishes in the integral, another change of variable gives 

f fdU" = cos{k^i) 

(37) /cosK/2)\ .cosK/2) 

^ / "''^ 1\ cos{k(p) ^f^dv? 

J-aj2 \ C0S(V3) y C0S2(V9) 

= cos{kipi)Akakizy). 

The secants Hz^,,pi, I — 1, . . . , n^^, intersect only at a finite number of points and 

(38) ^/(t) = X! cos(fc</7i)Aafc(z7)- 

Using Lemma [TT] to evaluate the sum, we arrive at the claim. □ 

Corollary 13. Let f be as in Lemma [71[ If E E and 'y E Te is such that 
i-y + K-y = 0, then 

(39) Gfi-f) = Skh)Akak{z^). 

In particular, if E = {0} so that i-y = = 0, for injectivity it is sufficient to 
consider only those 7 for which gcd{ny,m^) — 1. In this case 

(40) ^/(7) = %(-l)'^'"-/"- Aafc(^7) 
if k/riy is an integer and Gf{"/) = otherwise. 

Proof. The function / may be written as a component of a Fourier series in the 
variable d = "d — i^: 

(41) /(r, §) = hk{r) cos(y) + 6fc(r) sin(y). 

This corresponds to choosing coordinates so that t-^ — 0, so that the problem 
reduces to the situation in Lemma |12[ Due to elementary trigonometric identities 
the coefficients of this Fourier series relate to the ones in Lemma [5] so that 

dkir) = ak{r) cos(fct^) — bk{r) sin(fct-y) and 

(42) 

^fe(^) — afe(^) sin{kLj) + bk{r) cos(fct^). 
Lemma 1 1 2 1 indicates now that 

Gfil) = Ski-y^icosikK^^AkOkizj) + sm{kK^)Akbk{Z'y)) 
(43) = Ski^)icos{kLj)Akak{z^) - sm{kLj)Akbk{zj)) 

= Sk{j)Akak{zj), 
which is the first part of the claim. 



BROKEN RAY TOMOGRAPHY IN THE DISK 



11 



Let US now turn to the case E ~ {0}. Given 7, let 7 be such that = 
riy / gcd{n^ , m^) , m^y — / gcd{n^ , m^) , = a^, and tr(7) = tr(7). It is geo- 
metrically obvious that 7 contains gcd(n-y,m^) copies of 7 and so Qf{'^) — Qfi"f)- 
Thus it suffices to consider those 7 £ Tfo} with coprinie and m^. 

Now that n-ya-y = 27rrn^ and gcd(n^,r7i^) = 1, ka^ e 27rZ if and only if k/rij 
is an integer. Thus the definition of Sk{j) together with the result Gfij) — 
Skil)^k0.kiz'y) gives the second part of the claim. □ 

Note that the functions bk do not appear in the above formulas for Gfi'y)- This is 
due to the assumption that the broken ray is symmetric with respect to the normal 
to dD at angle 0, which causes functions antisymmetric with respect to this line to 
integrate to zero over the trajectory. 

4.2. Quasianalytic functions. The class of functions for which our argument 
shows the uniqueness result of Theorem [3] consists of functions uniformly quasian- 
alytic in the angular variable as defined below. There is, however, no implication 
that this is the optimal class in which the theorem should hold. 

Definition 14. Let M = (Af„)^Q be an increasing logarithmically convex se- 
quence of strictly positive real numbers. A sequence (flfe)^]^ of real numbers is 
defined to be in the class S*^ ( Af ) if there is a constant R > such that for each n 
it satisfies 

CO 

(44) ^fc'"|afe|' <M„i?". 

The class S'^{M) is called a quasianalytic class of sequences if the defining sequence 
M satisfies 



00 



(45) E 



: 00. 



ri=0 

A function / : M — > M with a period 2tt belongs to the class (M) if it can be 
written as a uniformly convergent Fourier series as 

CO 

(46) f{x) — E('^fc cos(fca;) + bk sm{kx)) 

k=0 

such that the sequences {ak)kLi {bk)kLi £^re in the class S'^{M). If the class 
S'^{M) is quasianalytic, then C^iM) is a quasianalytic class of functions. Func- 
tions in such a class are quasianalytic functions. 

If a continuous function / : I) — >■ R is written as a Fourier series as in Lemmajsj it 
is uniformly quasianalytic in the angular variable if for each R G (0, 1] the collection 
of functions {/(r, •) : r > R} belong to the same quasianalytic class. 

In dimensions higher than 2, we define a a function / : — > M, ?i > 3, to 
be uniformly quasianalytic in the angular variables if its restriction to any two 
dimensional plane intersecting the origin is uniformly quasianalytic in the sense 
defined in two dimensions. 

The class C"^(Af) of 27r-periodic smooth functions is a quasianalytic class of 
functions in the sense defined above if and only if it is so in the classical sense: if 
/ is in a quasianalytic class and the function and all of its derivatives vanish at 
any point, then the function is identically zero. For this result and more discussion 
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on periodic quasianalytic functions we refer to [6l Section V.2]. We only remark 
here that quasianalytic functions are smooth and all derivatives of the Fourier 
series (which can be taken term by term) converge uniformly to the corresponding 
derivative of the function. 

Angular differentiation of a function uniformly quasianalytic in the angular vari- 
able may also be carried out term by term and the series converges uniformly, 
provided that the angular derivatives of all orders satisfy the Dini-Lipschitz condi- 
tion. This is why we make the assumption in Theorem [3] 

The important property of quasianalytic functions that motivates the use of this 
class is the following lemma. 

Lemma 15. Let {xk)'^Q be a bounded sequence and z G (0,1). // a uniformly 
quasianalytic function f is written as a Fourier series as in Lemma^ and 

oo 

(47) J2^-ey'xkAkak{z) = 

fe=0 

for all n ^ 0,1, ... , then XkAkaki^) = for each k. 

Proof. The functions afc(r, •), r > z, belong to the same quasianalytic class C^(Af ). 



By Lemma 10 ;3) we have (Aafe(z))^o ^ S*{M) 



Define a function g : K — >■ E so that 

oo 

(48) g{x) = ^ XkAkUkiz) cos(fcx). 

fc=0 

Now g is quasianalytic and obviously its odd derivatives vanish at zero. By as- 
sumption also the function itself and its even derivatives vanish at zero, so g must 
be identically zero. This proves the lemma. □ 

5. Tomography from a singleton 

We prove Theorem [T] in the special case when the unknown function / satisfies 
the Dini-Lipschitz condition. This result is not needed for the proof of the general 



case (which was presented in Section 2.3), but we present it since it is simple and 



uses the tools developed for an open set of tomography. Both proofs are constructive 
and the reconstructions are based on the 'long trajectory limit', where n-y tends to 
infinity. 

Proposition 16. For f e C^^, Qf : Te E determines the zeroth Fourier 
component ao(r) = fdH^ when E is a singleton. 

Proof. As noted in Corollary [13] it is sufficient to consider those 7 G F^; for which 
gcd(7i^,m-y) — 1. By continuity of / it is enough to reconstruct it outside the 
boundary dD. 

Let Ofe and bk be as in Lemma [s] Let g{r, ??) = aQ{r) and 

K 

(49) fK{r, d) = y(ak{r) cos{M) + 6fe(r) sin(fc,?)). 



fc=i 



By Corollary 13 Ggij) ~ d ^.4oao(2:-Y), and using this we show that Gfi'j) — > Qg{l) 
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Let £ > 0. By Lemma|8]there is Kf, such that ||/ — 5 — /k, II < By Corollary 13 
QIkAi) = whenever 7 e F^^^^, whence 

(50) 1^/(7) - 5.9(7)1 = \G{9 + /a7 - /)(7) - GfKAl)\ < ^ + 0. 

Fix any w £ [0, 1). There is a sequence (7i)ieN such that n-y. > i, and z^. w 
(and thus d^. — !■ 2\/l — w'^) as i — > cx). The estimate obtained above gives thus 
linii_j.oo ^/5(7i) — 5(1 — w'^)~^^'^AQaQ{'w), allowing one to reconstruct the function 



^0^0 on [0, 1). Injectivity of the Abel transform (Lemma 10 1 proves the claim. □ 



By Corollary 13 and Theorem[l]we have the following for / e C^^ and a singleton 
E (chosen to be at angle zero): G/Ife — if only if oq = and 

(51) ^(-l)'='"-/"-Aafc(^7) = 

for all 7 G F^; with gcd(n.^,m^) — 1. Here the sum is over all k with n^\k. 

First, nothing can be said about the functions bk- Furthermore, there are a 
number of functions such that AkO-kiz-y) = whenever rt^jfc; by injectivity 
-4fe(Co([0, 1])) is an infinite dimensional subspace of C([0, 1]), and the space {/ G 
-4fe(Co([0, 1])) : f{zj) = when nj\k} has finite codimension in Ak{Co{[0, 1])) (for 
each k, there are only finitely many 7 with nj\k) and therefore infinite dimension. 
Hence the condition G/Ite = leaves great freedom to the functions Ofc, k > 1, 
and we can efficiently recover only oq from G/Ife- 

One can only possibly image the unknown function / in the set IJ^gp^ ^1(7). 
When E is countable, this set is dense in the disk but has Hausdorff dimension one, 
so without continuity assumptions no recovery results can be expected. 

6. Tomography from an open set 

We denote by (7"'™ = C'^jjC™g the space of real valued functions on the punc- 
tured disk D* — D \ {0} which have n continuous partial derivatives with respect 
to the radial and m to the angular variable. The angular derivative is denoted by 
dang- We require no regularity at the origin and correspondingly restrict Gf to 
= {1 ^ : ^ tr(7)} for / G C"'™. Including or excluding the origin in the 
domain of the unknown function does not alter the results of Theorems [T| and [3] 
besides reconstructibility at the origin. 

Let a rotation operator G M, act on functions so that y\^f{r, •&) = f{r,'d + 

if). Let another rotation operator act on trajectories instead of functions so 
that L.yi^'Y = i~f + ip, Km.^-y — Kry + ip and other parameters describing 7 are unaltered. 
Geometrically Ui^p simply rotates the trajectory by an angle (p. 

We can only show the uniqueness result of Theorem [3] for functions uniformly 
quasianalytic in the angular variable, but if it would hold for functions in (7"'°°, 
the following lemma implies that it would also hold for the class C"'*'. 

Lemma 17. For E C dD open and any n, the following hold true: 

(1) If the uniqueness result of Theorem^would hold for the class C"'°°, then 
it would also hold for C"'". 

(2) If f e C"^" and gf{j) = for all 7 G F|j, then also Gdl^J{j) = for all 
7 G T*^ and i G N, z < m. 

(3) Suppose f G C"^" and gf{-f) = for all 7 G F^j. Fix j e T*^ and S > 
such that yiip"f G F^ whenever \ip\ < S. Then G^tpfi'^) = for all \ip\ < S. 
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Proof. Each part is based on the following observation. The change of the integral 
of / over a trajectory is the same if the function is rotated clockwise as if the 
trajectory is rotated counterclockwise by the same angle ip: 

(52) gf{^-^i)^g^^fh). 

(1) Let us denote a shrinked set of tomography by 

(53) = {a e E : d{x,dD\E) > S}. 

For a sufficiently small S > the set E^ is nonempty and open. 

Let ?7 : M M be the standard moUifier and let — e~^ri{x/e). Let e < S. 

If f e C"'" is such that g/lrj, = 0, then gf{m^-/) = for all 7 e T*^, and 
\ip\ < S. Thus for any 7 G T*^^ we have 

/>oo 

0- / %(^)e/($H^7)d^ 

(54) 



— 00 
00 



7j,i^)gm.^f{j)dip 

where the convolution is with respect to the angular variable of /. Therefore 
gf\r'^ = implies that g{r], * /)|r^, - 0. If / G C"'0, then r,, * f e C"^°°. 

If Theorem jsj holds in the class C"'°°, then 5(?7e * /)|r'j =0 implies that 
rie * f = 0- Since this holds for all e e (0, (5) and i]e * f ^ f uniformly as e 0, we 
find that / = 0. 

(2) Suppose TO > 1 and / £ C"'"'. It is enough to show the result for the first 
derivative, from which the argument can be iterated. 

Dividing Eq. ( 52 1 by and passing to the limit — > we obtain 

(55) 9^5/(9^^7)1^.0 = -5d.nJh) 

for all 7 G r*^. Since gf vanishes on T*^, this implies that gdangf does, too. 

(3) Trivial. □ 

Proof of Theorem]^ We show that ^7(7) = for all 7 e Lb impHes / = if / 
is uniformly quasianalytic in the angular variable and satisfies the Dini-Lipschitz 
condition. 

Let / be written as a Fourier series as in Lemma [8] By Theorem [T] we already 
know that ag = 0. Differentiating term by term, we have 



00 

(56) dllj{r,d) = Y,{-er{ak{r)cos{M) + hk{r)sm{kd)) 

k=0 



for any n £ N. We may choose coordinates so that the angle zero lies in E. Fix 
any 7 g F^; such that + = 0; such a trajectory can be found since E is open. 

By Corollary |13[ Lemma 17 and the assumption of vanishing broken ray trans- 
form we have 



00 

(57) ^(-fc2)"5,(^)^,afe(z^) = 

k=0 
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for n = 0, 1, . . . . Hence by Lemmas 11 and 15 we have that Sk{j)Akak{z-y) = oj^ 

If = —K-y ^ ttQ. we have Skij) 7^ for all k by Lemma 11 and so Akikizj) — 
for all k. The set of such z-y that = —Kj ^ ttQ and 7 S F^; is dense on [0, 1], so 
continuity of AkCik (Lemma |10[ ) implies that AkCik = 0. This, by the injectivity of 
Ak (Lemma 10 1, implies that each is identically zero. 

We have thus demonstrated that the function / is antisymmetric with respect to 
the normal to dD at angle zero. The same conclusion may be drawn with any line 
intersecting E and the origin. If / is antisymmetric with respect to two different 
lines going through the origin, / vanishes everywhere. □ 

Notice that by Proposition [s] / G C with Qf = vanishes in f2 \ co{dil \ E). 
If / is also quasianalytic in the angular variable, it follows that / vanishes outside 
the open disk of radius (i(0, Cl \ co{dVl \ E)) centered at the origin. The fact that / 
also vanishes inside this disk, as shown above, is nontrivial. We also remark that 
we have not employed the support theorem in our proof of Theorem |3] 

Proof of Corollary^ Take any two dimensional subspace F of M" such that it 
meets E. As noted in the proof of Corollary |2] at the end of Section |2.3[ one ends 
up with the two dimensional broken ray transform in the disk n F. Since EOF 
is open in the circle S"^^ F, Theorem [s] guarantees that Q f restricted to those 7 
that lie in F uniquely determine / in B" n F. Since M" can be written as a union 
of such subspaces F, Qf determines / everywhere. □ 
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This conclusion may also be drawn as follows. Write / as a complex Fourier series: /(r, 1?) = 
Efcez Define h : R ^ C by h{v^) = 5^^/(7) = Efesz (7)^|fe| ^fe (^7)- The 

sequences {Sk{'y)A.kO.±k{r))'^Q are quasianalytic sequences, whence h is a quasianalytic function. 
By part (3) of Lemma |l7| fe vanishes on some open interval, and so by quasianalyticity it vanishes 
everywhere. Thus each Fourier coefficient (7)^1^,1 aj. (2:-^) is zero. 



